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Darwin (1920) noted that when radiation can be neglected it should be possible to eliminate 
the radiation degrees-of-freedom from the action of classical electrodynamics and keep the discrete 
particle degrees-of-freedom only. Darwin derived his well known Lagrangian by series expansion in 
v/c keeping terms up to order (v/c) 2 . Since radiation is due to acceleration the assumption of low 
speed should not be necessary. A Lagrangian is suggested that neglects radiation without assuming 
low speed. It cures deficiencies of the Darwin Lagrangian in the ultra-relativistic regime. 
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When radiation can be neglected the Lagrangian of 
classical electrodynamics, putting (3 — v/c, can be writ- 
ten, 



L = J2\ -™«c 2 yl -Pl + Y \Pa ■ Mr a) - <t>{r a )] 
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In 1920 Darwin [1] expanded the Lienard-Wiechert po- 
tentials to second order in /3 — v/c and thus found that, 
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and (hats are used for unit vectors), 
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give the correct Lagrangian to this order. More recent 
derivations can be found in a few textbooks [2-4] . In par- 
ticular Jackson [4] notes that using the Coulomb gauge 
(V • A = 0) makes the electrostatic Coulomb potential 
4> exact and moves all approximation to the vector po- 
tential A which obeys the inhomogeneous wave equation 
with the transverse (divergence free) current as source. 
The Darwin Lagrangian results when the term c~ 2 d 2 /dt 2 
in the d'Alembert operator is neglected so that the equa- 
tion becomes a Poisson equation. 

The Darwin Lagrangian has been shown to be a 
unique approximately relativistic Lagrangian (Woodcock 
and Havas [5], Kennedy [6]). It can be derived from 
the Fokker-Wheeler-Feynman action-at-a-distance the- 
ory (Anderson and Schiminovich [7]), and it is useful in 
various fundamental studies of electrodynamics [8-11]. 
The magnetic interaction described by the Darwin La- 
grangian is essential in relativistic many-electron calcu- 
lations as noted by Breit and others [12-15]. It has found 
applications in nuclear physics [16, 17], and especially in 
plasma physics, for numerical simulation [18-22], ther- 
modynamics and kinetics [23-27], as well fundamental 
theory [28-30]. Barcons and Lapiedra [31] noted that 
the Darwin approach is not valid for a relativistic plasma 



and therefore used a different approach to its statistical 
mechanics. 

Corrections to the Darwin Lagrangian have been dis- 
cussed. Since a system of particles with identical charge 
to mass ratio does not dipole radiate a higher order ex- 
pansion should be meaningful for such systems [32-34]. 
To that order, however, acceleration inevitably enters 
and must be handled in some way. Others have argued 
that since radiation is due to acceleration, v/c expansion 
is irrelevant, and further that radiation can be negligible 
even if the particle speeds are considerable (Trubnikov 
and Kosachev [35], Frejlak [36]). We will pursue that 
lead here. 

One frequently encounters the statement that the Dar- 
win approach neglects retardation. This may be due to 
the fact that the, nowadays best known, elegant deriva- 
tion by Jackson [4] hides the complications due to re- 
tardation. Nevertheless it is wrong. The derivations by 
Darwin [1] and by Landau and Lifshitz [2] show that the 
contribution of retardation to the Coulomb potential in 
the Lorenz gauge, is quite large. The main acceleration 
dependent part, however, vanishes either, as in Darwin's 
derivation, because it gives a total time derivative term 
in the Lagrangian, or, as in Landau and Lifshitz, because 
of a gauge transformation (to the Coulomb gauge) . Both 
these derivations also show that the velocity dependent 
part of the retardation is handled exactly to order (v/c) 2 . 

A natural idea that does not work is to assume constant 
velocities and use the known exact Lienard-Wiechert po- 
tentials for that case in (1). Darwin's original deriva- 
tion shows that this does not give the electric interac- 
tion to sufficiently accuracy. It is important to note that 
gauge invariance (for a review, see Jackson and Okun 
[37]), which is valid for the exact theory, does not nec- 
essarily hold for approximations. We therefore impose 
the Coulomb gauge (for a recent discussion see Heras 
[38]) and then solve the inhomogeneous wave equation for 
A assuming constant velocities in the transverse current 
density. In this way one treats the electric interaction 
exactly, neglects acceleration in the solution for A, but 
do not assume low speeds. 
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The constant velocity exact Coulomb gauge vector po- 
tential does not seem to be well known. A special case 
was solved by Labarthe [39]. The explicit general solu- 
tion has recently been published by Hnizdo [40] who used 
a gauge transformation function given by Jackson [41] to 
find it, starting from the corresponding Lienard-Wiechert 
potentials. Denote by, R = r — r'(t), the vector from the 
source particle at r'(t), with charge e, to the field point 
r, so that (3 = r\t)/c. If we then put, 



r) = R x (3, 



(4) 



Hnizdo's solution, which assumes the source particle to 
be at the origin at time t — and to have constant ve- 
locity along the x-axis, i.e. r'(t) — cfltx, can be written, 
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P<f>L - (4>L ~ <t>c)/P, 
y 2 + z 2 ( ( l } L-^c)/0, 
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at t = 0, so that x,y,z are the components of R. 
Here, (f>c — e/R, is the Coulomb potential, and, (j>L — 
<Ac/-\/l — ? 7 2 > Lorenz gauge form. One notes the iden- 
tity l/(y 2 + z 2 ) = {l3/rj) 2 /R 2 . Using this, and that the 
only relevant vectors are R and /3, one can, by expressing 
everything in terms of these, or scalar and vector prod- 
ucts involving these, arrive at the coordinate independent 
form, 

e \g{rf)l3 + h{rf){l3-R)R 
Ac(r) = -i ± (8) 



R 

of Hnizdo's solution. Here we have introduced the nota- 
tion, and the functions g and h are defined by, 
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It is remarkable that an equivalent vector potential has 
been found by Crater and Lusanna [9] in a canonical for- 
malism. When the momenta (denoted k) of Eq. (5.28) 
of [9] are replaced by m/3/y/l — ft 2 the expression (8) 
is recovered. The authors of [9] use a relativistic phase 
space formalism and assume that charges are anticom- 
muting Grassmann variables. In this way they treat the 
Pauli exclusion principle semiclassically. Their Hamilto- 
nian formalism, which must entail the neglect of acceler- 
ation in an indirect way, is mainly intended for treatment 
of bound states. The Hamiltonian based on the ordinary 
Darwin Lagrangian (1) is discussed in [25]. 

The explicit expression for the interaction Lagrangian 
of two particles that results when (8) replaces (3) in (1) 



is, 
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where, rj 2 = (r a & x /3 a ) 2 . We now consider two special 
cases. If the velocity of a particle is parallel to the inter- 
particle vector to another particle, rj = 0, so the Darwin 
interaction needs no correction in these cases. Assuming 
that two particles have equal velocities Vi = v 2 = v 
parallel to f 2 i we find the the interaction term in (1) 
gives, 
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and, 



and (11) gives the same result. One sees that this term, 
and the corresponding force, goes to zero in the ultra- 
relativistic limit v — > c. Now consider instead the inter- 
action of two particles that move with equal velocities, 
Vi = v 2 = v, but side by side, so that, r 2 i -L v. The 
interaction part of the Lagrangian (1) is then, 



h(x) 
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(10) 



Note that g(l) = 1 but that h diverges for x = 1. From 
these expansions it is clear that the leading terms give 
the vector potentials (3) of the Darwin Lagrangian. 

The vector potential of the original Darwin Lagrangian 
is thus recovered from (8) when r] = 0. One notes that 
in the derivation of (3) there was no need to assume that 
the velocity is constant since the solution to the Poisson 
equation does not require retardation, while it is neces- 
sary for solving the wave equation. Hence the assumption 
of constant velocity. One purpose of a Lagrangian is, af- 
ter all, to find equations of motion that determine the 
accelerations. If it is necessary to know them beforehand 
a Lagrangian approach is pointless. 



1 e x e 2 v 2 e x e 2 
^12 — 77- 
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One sees that even in the limit v — > c the Coulomb in- 
teraction dominates; the magnetic interaction can only 
compensate for half of it. This is clearly wrong, however. 
It is well known that in an ultra-relativistic beam the 
transverse Lorentz force cancels the transverse Coulomb 
repulsion (see e.g. [42]). 

Let us instead use the vector potential (8) and the 
corresponding interaction (11). We first note that in this, 
side by side, case r/ 2 = v 2 jc 2 and thus r/ 2 = 1 in the limit 
v — > c. In this limit g(l) = 1 and h diverges. The two 
scalar products in the second term will, however, be zero 
and a simple investigation shows that this compensates 
for the divergence of h, so that term does not contribute. 
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Finally we get, 



r 111 2N e l e 2 v 
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(14) 



and in the limit v — > c this term and the corresponding 
force, are zero, as they should, when (11) is used. 

In conclusion, the Lagrangian obtained by using the 
exact constant velocity Coulomb gauge vector potential 
(8), instead of the A used in (3), has been derived with- 
out assuming that v/c is small, only that accelerations 
are not needed in estimating the Coulomb gauge vector 
potential. In this way all velocity dependent retardation, 
and, as discussed above, also the main part of the ac- 
celeration dependent retardation, is accounted for. We 
have also shown that using this Lagrangian we account 
correctly for the pinching of an ultra-relativistic beam, 
something the original Darwin Lagrangian does not do. 
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